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Abstract 


A  combined  experimental  and  modeling  study  of  neat  and  NH3 -doped  (^ = 1),  30-Toit  flames 
is  reported.  The  major  species  concentrations  are  measured  by  molecular  beam  mass 
spectrometry  (MB/MS),  whereas  the  minor  species  OH,  NH,  and  0-atom  concentrations  are 
measured  by  laser-induced  fluorescence  (LIF).  The  species  NO  is  measured  both  by  LIF  and 
MB/MS,  and  Oj  by  MB/MS.  The  flame  temperatures  are  measured  both  by  OH  and  NH  LIF  and 
by  thin  wire-thermometry.  The  flames  are  modeled  with  PREMIX  using  the  temperature  profiles 
and  several  detailed  chemical  mechanisms  as  input.  The  mechanisms  include  the  GRI  2.11, 
SSLA,  and  their  derivatives.  The  SSLA  mechanism  was  developed  previously  in  our  laboratory 
from  a  critical  literature  review.  Calculations  using  all  the  mechanisms  predict  fairly  well  the 
profiles  of  the  major  species  for  both  neat  and  doped  flames.  However,  both  the  SSLA  and  GRI 
2.11  calculations  fail  to  predict  the  postflame  Oj  concentration  in  the  neat  flame,  the  drop  in  the 
O2  concentration  with  the  addition  of  NH3,  and  the  NH3  decay  in  the  doped  flame.  Sensitivity 
analyses  suggest  refinements  to  the  SSLA  and  GRI  2.11  mechanisms.  The  experimental  results 
are  predicted  rather  well  using  a  modified  SSLA  mechanism  in  which  the  NH+N0=N20+H 
reaction  rate  is  decreased  and  the  N20+M=N2+0+M  reaction  rate  and/or  H2O  third  body 
efficiency  is  increased  to  the  limit  of  their  uncertainty.  Rate  analyses  performed  on  the  modeled 
calculations  reveal  the  reactions  important  to  NO,  O2,  NH,  OH,  and  O-atom  production  and 
consumption  and  NH3  consumption.  These  reactions  are  presented  and  discussed. 


u 


Acknowledgments 


We  thank  Drs.  W.  Anderson,  A.  Kotlar,  and  R.  Pastel  of  the  U.S.  Army  Research  Laboratory 
(ARL),  Dr.  D.  Dayton  of  the  National  Renewable  Energy  Laboratory,  and  Prof.  G.  Singh  of  the 
University  of  Maryland  Eastern  Shore  for  their  helpful  discussions.  We  also  thank  Dr.  W.  Anderson 
for  carefully  reviewing  this  manuscript.  This  work  was  supported  by  the  ARL  Mission  Program  on 
Combustion  and  the  Strategic  Environmental  Research  and  Development  Program  (ERDP)  on 
cleanup  (ARMY-7 1 3-94  and  ARMY-729-94).  Support  from  the  National  Research  Council  (NRC) 
ARL  Post-Doctoral  Fellowship  Program  (D.  T.  Venizelos),  and  the  Productivity  Capital  Investment 
Program  (R.  C.  Sausa)  is  gratefully  acknowledged. 


m 


Table  of  Contents 


Acknowledgments  . 

List  of  Figures . 

List  of  Tables  . 

1.  Introduction . 

2.  Experimental  . 

3.  Modeling . 

4.  Results  and  Discussion . 

5.  Conclusion  . 

6.  References . 

Distribution  List . 

Report  Documentation  Page 


Paee 

iii 

vii 

ix 

1 

2 

6 

8 

25 

27 

29 

31 


V 


List  of  Figures 


Figure  Cage 

1 .  Experimental  (...)  and  Fitted  (-)  NH  Spectra  in  the  302.2-302.7-nm  Region  at 

5  mm  for  a  30-Torr  Hj/NjO/Ar  Hame  Doped  With  NH3  .  4 

2.  Temperature  Profiles  for  Both  Neat  and  NHj-Doped  Flames .  6 

3.  Plot  of  NH3  Sensitivity  Coefficients  as  a  Function  of  Height  Above  Burner  for 
Various  Reactions  Using  Both  SSLA  (Top)  and  GRI  2.11  (Bottom) 

Mechanisms . 

4.  Experimental  and  Modeled  NH3  Profiles .  15 

5.  Experimental  and  Modeled  O2  Profiles .  16 

6.  O2  Reaction  Rates  From  SSLA  Modified  Calculations  for  Both  Neat  and 

NH3-Doped  H2/N20/Ar  Flames .  17 

7.  Experimental  and  Modeled  NO  Profiles  for  Neat  and  NH3-Doped  H2/N2O/ Ar 

Flames  . 

8.  Reaction  Rates  of  NO  From  SSLA  Modified  Calculations  for  Both  Neat  and 

NH3-DopedH2/N20/Ar  names .  20 

9.  NH  Experimental  and  Modeled  Profiles  for  Neat  and  NH3-Doped 

names  .  ^1 

10.  NH  Reaction  Rates  From  SSLA  Modified  Calculations  for  Both  Neat  and 

NH3-DopedH2/N20/Ar  names .  23 

1 1 .  OH  Experimental  and  Modeled  Profiles  for  Neat  and  NH3-Doped  H2/N20/Ar 

names  .  24 

12.  0-Atom  Experimental  and  Modeled  Profiles  for  Neat  and  NH3-Doped  H2/N20/Ar 

names  . 


vu 


List  of  Tables 


Table 


Page 


1.  Sensitivity  Coefficients  for  O2  in  a  Neat  and  NH3-Doped  Flame  at  15.0  mm 

Above  the  Burner  Surface . 

2.  Sensitivity  Coefficients  for  NH  in  a  Neat  and  NH3-Doped  Flame  at  4.1  and 

5.3  mm.  Respectively,  Above  the  Burner  Surface . 


ix 


1.  Introduction 


The  H2/N2O  flame  system  has  received  considerable  attention  in  recent  years  because  it  is  an 
appropriate  analog  for  experimental  and  modeling  efforts  to  study  NO^  pollutant  formation  and 
nitranninp.  propellant  combustion  [1—3].  Such  a  flame  system  can  provide  additional  information  on 
the  fundamental  reaction  mechanisms  relevant  to  the  nitrogen  chemistry  in  more  conc5)lex  combustion 
systems.  Recently,  the  role  of  NH3  in  converting  NO  to  final  products  was  investigated  in  our 
laboratory  using  neat  and  NH3-doped  H2/N20/Ar  low-pressure  flames  (O'- 1.1)  [4].  The  profiles  of 
several  species,  such  as  Hj,  N2O,  NH3,  HjO,  Nj,  and  NO,  in  both  neat  and  doped  flames,  were 
obtained  with  molecular  beam  mass  spectrometry  (MB/MS).  The  experimental  profiles  were 
compared  to  predicted  profiles,  generated  with  the  one-dimensional,  laminar  flame,  PREMIX 
computer  code  using  a  detailed  chemical  mechanism  ( SSLA)  derived  from  a  critical  literature  review 
[4].  The  modeled  profiles  of  the  major  species  agreed  very  well  with  the  experimental  profiles  for 
the  neat  flame  and  reasonably  well  for  the  doped  flame.  However,  there  were  discrepancies  between 
the  experimental  and  predicted  NO  profiles.  Specifically,  the  model  predicted  a  decay  of  the  NO 
profile  in  the  postflame  region  for  the  doped  flame,  while  the  experimental  profiile  exhibited  a  plateau. 
In  addition,  the  model  overpredicted  the  NH3  mole  firaction  in  the  postflame  region  and  that  this 
overprediction  was  the  primary  cause  of  the  predicted  NO  decay.  Calculations  indicated  that 
refinements  to  the  chemical  mechanism  used  by  the  numerical  model  were  necessary. 

Such  refinements  require  comparing  not  only  the  experimental  and  predicted  profiles  of  major 
species  but  also  con^aring  the  profiles  of  intermediate  ^ecies.  Unfortunately,  these  measurements 
prove  to  be  difficult  with  MB/MS  for  a  H2/N2O/NH3  system  because  of  interference  between  isobaric 
species  such  as  O  and  NH2,  or  OH  and  NH3.  Also,  the  measurements  of  NH  may  be  obscured  by  the 
interference  from  NH2  and  NH3  because  of  ionization  fi^gmaitation  effects. 

Reported  in  this  paper  is  a  combined  experimental  and  modeling  study  of  neat  and  NH3-doped 
H2/N2O  flames.  The  current  work  extends  the  previous  study  in  our  laboratory  by  introducing 
additional  measurements  of  the  O,  OH,  NH,  and  NO  species  using  laser-induced  fluorescence  (UP) 
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and  Oj  using  MB/MS .  Rame  temperatures  are  measured  both  by  thin-wire  thermometry  and  OH  and 
NH  T.TF.  The  flames  are  modeled  using  the  PREMIX  flame  code  with  the  detailed  SSLA  [4]  and 
GRI  2.11  [5]  mechanisms  and  their  derivatives  as  input.  Rate  and  sensitivity  analyses  reveal  key 
reactions  in  the  mechanisms,  as  well  as  the  pathways  for  species  production  and  consumption. 


2.  Experimental 

The  experimental  apparatus  consists  of  a  variable-pressure  burner  equipped  for  MB/MS,  LIF, 
andthin-wire  thermometry.  The  details  are  discussed  elsewhere  [1,2, 4].  Briefly,  both  the  neat  and 
NHj-doped  H/NjO/Ar  flames  were  supported  on  a  6-cm,  stainless  steel,  McKenna  flat  burner 
mounted  in  a  cylindrical  vacuum  chamber  maintained  at  30  Torr.  The  neat  flame  was  stabilized  by 
flowing  a  mixture  of  Hj,  NjO,  and  Ar  at  1 .41 , 1 .43,  and  1 .06  slm,  respectively.  This  flame  was  then 
doped  with  0. 1 4  slm  NH3.  The  experimental  apparatus  was  reconfigured  for  the  present  smdy  so  that 
the  T.TF  and  thermocouple  flame  measurements  were  performed  at  the  same  point  and  at  the  center 
of  the  burner  directly  beneath  the  MB/MS  probe.  A  comparison  between  the  MB/MS  and  UF 
concentration  profiles  of  NO  indicated  that  the  MB/MS  sampling  takes  place  approximately  1 .5  mm 
(five  nozzle  diameters)  below  the  probe  tip.  Cattolica,  Yoon,  and  Knuth  [6]  also  observed  a  shift 
of  five  nozzle  diameters  when  measuring  the  OH  radical  by  both  laser  absorption  spectroscopy  and 
MB/MS  in  a  stoichiometric  CHyair  flame  [6].  The  present  experimental  configuration  allows  for 
a  more  accurate  measurement  of  the  flame  temperatures  and  species  concentrations  because  the  LEF 
and  thermocouple  measurements  are  recorded  exactly  where  the  MB/MS  is  sampling.  All  the 
MB/MS  data  are  shifted  ~1.5  mm  toward  the  burner  surface  for  direct  comparison  to  the  UF  and 
thermocouple  data. 

A  20-Hz,  Nd-YAG  (Quantel,  YG581C)  pumped  dye  laser  (Lambda  Physik,  FL3002)  witha  BBO 
frequency-doubling  crystal  was  used  for  the  NO,  OH,  NH,  and  0-atom  IF  measurements.  Pulse 
energies  of  100-500  pJ  with  5%  shot-to-shot  fluctuation  and  10-ns  duration  were  obtained  after 
frequency  doubling.  The  2-mm  diameter  laser  beam  was  focused  with  a  200-mm  focal  length  lens 
for  0-atom  UF  and  a  500-mm  focal  length  lens  for  NO,  OH,  and  NH  UF.  A  shorter  focal  length 
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lens  was  required  for  0-atom  UF  because  of  the  multiphoton  nature  of  the  excitation  process.  A  plot 
of  LEF  signal  as  a  function  of  laser  energy  showed  a  near  quadratic  dependence,  suggesting 
two-photon  excitation. 

* 

The  excitation/detection  schemes  for  profiling  the  aforementioned  species  involve  die  following: 
(1)  excitation  of  NO  (0,0)  Qj  (26.5)  and  Rj  (25.5)  transitions  at  225.37  [7]  and 

225.03  nm,  respectively,  with  subsequent  (0,1)  emission  detection  at  236  nm;  (2)  two-photon 
excitation  of  the  O  3p^P  *-  2p^P  transition  near  226  nm  with  3p^P  -  3s^S  emission  detection  at 
845  nm;  (3)  excitation  of  OH  A^S"  -  X^II  (1,0)  Ri  (7.5)  transition  at  281.1  nm  with  (1,1)  emission 
detection  at  313  nm,  and  (4)  excitation  of  the  NH  A^S*  -  X^H  (1,0)  R2  (9)  transition  at  302.65  nm 
with  (1,1)  emission  detection  at  337  nm.  The  fluorescence  was  collected  at  an  angle  normal  to  the 
excitation  laser  beam  and  viewed  with  a  photomultiplier  tube  (Hamamatsu  R955)  equipped  with  an 
appropriate  10-nm  (fwhm)  interference  filter  for  O  and  OH  detection  or  a  1/8-m  wide-slit  (~10  nm 
at  fwhm)  monochromator  (Jarrell  Ash)  for  NO  and  NH  detection.  The  excitation  transitions  were 
carefully  selected  for  their  temperature  insensitivity  by  using  the  following  equation  [8]: 


J2  +  J  - 


hcB 


v; 


V  =  0, 


(1) 


where  J  is  the  rotational  energy  level,  T^^  is  the  average  flame  temperature,  and  k,  h,  c,  and  the 
usual  spectroscopic  constants.  Boltzmann  calculations  show  that  for  a  5%  change  in  a 2,000  Kflame 
temperature  the  UF  signal  changes  <1%.  All  of  the  species  profiles  were  obtained  with  a  10-ns 
gated  integrator  and  1,000-shot  averaging  to  minimize  the  errors  due  to  quenching  effects  and  laser 
intensity  fluctuations.  The  sampling  gate  was  positioned  near  the  signal  peak. 

The  flame  temperatures  were  measured  by  NH  and  OH  UF.  The  spectra  were  recorded  in  a 
linear  energy  regime  with  a  10-ns  gate  width  and  fitted  using  a  multiparameter  fitting  computer 
program  based  on  a  Boltzmann  rotational  population  distribution  analysis  [9, 10].  Figure  1  shows 
a  typical  UF  excitation  spectrum  of  the  R  branch  of  the  NH  A^S”"  -  X^II  (1,0)  transition  recorded 
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Figure  1.  Experimental  (...)  and  Fitted  (-)  NH  Spectra  in  the  302.2-302.7-nm  Region  at  5  nun 
for  a  30-Torr  Hj/NjO/Ar  Flame  Doped  With  NH3. 


for  the  NHj-doped  Hj/NjO/Ar  flame.  The  best  fit  to  the  NH  UF  spectram  is  shown  with  a  solid  line 
and  yields  a  temperature  of  1,867  K  ±40  K  (2SD)  at  5  mm. 


The  flame  temperatures  were  also  measured  with  a  thin  Pt-Rh  (6%)/Pt-Rh  (30%)  thermocouple 
that  was  coated  with  a  beryllium  oxide  ( 1 5%)/yttrium  oxide  mixture  to  avoid  surface  catalytic  effects 
[11].  They  were  corrected  for  radiative  losses  using 


T  -  T 

flame  meas 


e<^|(T^eas 


-  Tj). 


(2) 


where  and  T„  are,  respectively,  the  flame,  thermocouple,  and  ambient  temperatures,  e  is 

the  diermocouple  emissivity,  a  the  Stefan-Boltzmann  constant,  and  h  the  convection  heat-transfer 
coefficient  between  the  spherical  bead  of  the  thermocouple  junction  and  the  flame  gases.  The 
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convection  coefficient  is  calculated  using  the  flame  temperatures,  composition  of  the  flame,  and  gas 
velocities  for  different  heights  above  the  burner.  Calculations  reveal  that  the  correction  for  the 
radiative  losses  is  -350  K  in  the  postflame  region,  with  an  experimental  uncertainty  of  ~1 50  K.  As 
a  result  of  this  relatively  high  uncertainty,  the  ratio  eo/h  was  determined  for  several  different  heights 
above  the  burner  surface  using  the  rotational  temperatures  obtained  from  NH  UF  measurements  in 
the  region  near  the  burner  and  the  OH  UF  measurements  in  the  postflame  region  and  then  applied 
to  correct  the  thermocouple  temperature  profiles  in  both  the  neat  and  NH3-doped  flames.  The 
corrected  thermocouple  temperature  profiles  were  then  fitted  with  a  sigmoid-type  function  and  used 
in  the  PREMIX  calculations.  The  corrected  thermocouple  temperature  profiles  for  the  neat  and 
NHj-doped  flames,  as  well  as  the  temperatures  obtained  with  the  UF  technique  are  shown  in 
Figure  2.  The  overall  uncertainty  in  the  temperature  measurements  is  ±2  and  5%  (2SD)  for  NH  and 
OH  UF,  respectively.  As  shown  in  Figure  2,  when  3.5%  of  the  total  flow  of  NH3  is  added  to  the  a 
neat  Hj/NjO/Ar  flame,  the  temperature  profile  is  shifted  approximately  1  mm  away  from  the  burner 
surface  compared  to  that  of  the  neat  flame.  Also,  the  maximum  decrease  in  flame  temperature  with 
the  NH3  additive  is  -100  K  and  occurs  near  the  burner  surface,  while  there  is  no  change  in  the 
postflame  temperatures.  NAS A-Lewis  equilibrium  calculations  show  that  when  NH3  is  added  to  the 
neat  flame,  the  adiabatic  flame  temperature  decreases  by  19  K  [12]. 

The  species  profiles  for  both  neat  and  NH3-doped  flames  were  obtained  with  an  in-line  triple 
quadrupole  MB/MS  (ABB  Extrel  C50  TQMS)  operating  at  ionization  energies  from  15-30  eV, 
depending  on  the  species  of  interest,  to  minimize  interference  firom  other  species.  The  electron 
emission  current  was  maintained  at  0.1  ±0.01  mA  for  all  the  measurements.  The  profile 
measurements  of  the  stable  species  were  quantified  by  introducing  calibrated  mixtures  of  the  species 
of  interest  and  Ar  in  the  burner  chamber  and  measuring  the  ratio  of  the  species  to  Ar  signal 
intensities  at  room  temperature.  With  the  mass  spectrometer  settings  kept  constant,  flie  same  ratio 
of  intensities  was  tihen  measured  in  the  flame.  The  ratio  of  the  species  mole  fraction  to  the  Ar  mole 
fraction  in  the  flame  was  obtained  using  the  following  equation: 


fx.] 

fx,] 

'0 

/  \ 
Iat 

flame 

RT 

flame 

5 


2200 


HEIGHT  ABOVE  BURNER  (mm) 

Figure  2.  Temperature  Profiles  for  Both  Neat  and  NHj-Doped  Flames.  The  Flame 
Temperatures  Are  Measured  by  OH  LBF  for  the  Neat  Flame  (•),  NH  LIF  for  the 
Doped  Flame  (o),  and  Thin-Wire  Thermocouples  for  Both  Neat  (-)  and  Doped  (~) 
Flames.  The  Thermocouple  Measured  Temperatures  Are  Corrected  for  Radiation 
Losses  Using  the  LIF  Temperatures. 


where  X;  and  1;  are,  respectively,  the  mole  fraction  and  the  detected  signal  intensity  for  the  species 
of  interest.  The  signal  intensity  is  related  to  the  species  mole  fraction  by  =  Kj  Xj,  where  Kj  is  an 
instrument-dependent  sensitivity  factor.  Equation  (3)  is  derived  by  assuming  that  (Kj/Kat)  at  room 
temperature  (RT)  equals  at  the  flame  temperature  (flame).  The  validity  of  equation  (3)  was 

verified  in  previous  work  using  calibrated  mixtures  of  He  and  Ar  measured  under  ambient 
temperature  and  flame  conditions  [1, 2, 13].  Equation  (3)  neglects  Mach  number  focusing  that  is 
only  significant  for  Hj.  The  uncertainty  in  the  MB/MS  measurements  is  10%. 


3.  Modeling 


The  flames  were  modeled  using  the  Sandia  National  Laboratories  (SNL)  flame  code  PREMIX 
(ver.  2.55)  [14]  that  employs  the  CHEMKIN-II  (ver.  4.9)  chemical  kinetics  and  multicomponent 
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transport  libraries  [15,16],  with  the  measured  temperature  values  and  flow  rates  as  input  to  the  flame 
code,  and  with  several  chemical  mechanisms  involving  N-O-H  reactions.  Two  mechanisms  served 
as  the  basis  for  all  the  mechanisms  used  in  the  calculations.  The  first,  which  is  denoted  as  SSLA, 
was  developed  in  our  laboratory  and  was  published  in  Sausa  et  al.  [4] .  It  consists  of  87  reactions  and 
20  species  with  rate  constants  obtained  from  a  critical  literature  review.  The  second  mechanism  used 
is  a  subset  of  the  benchmark  and  popular  GRI  2.11  mechanism,  which  was  developed  for  the  Gas 
Research  Institute  for  natural  gas  ignition  and  flame.  The  GRI  2.11  mechamsm  contains  the  C-H-0 
chemistry  and  it  includes  the  nitrogen  chemistry  relevant  to  natural  gas  chemistry  and  rebuming.  All 
of  the  carbon-containing  species  and  associated  reactions  were  removed  from  this  mechanism, 
yielding  a  subset  mechanism  consisting  of  69  reactions  involving  19  species.  The  two  striking 
differences  between  the  SSLA  and  GRI  2.11  mechanisms  are  that  (1)  the  SSLA  does  not  contain  the 
NO2  and  H2O2  species  that  are  present  in  the  GRI  2.11  mechanism  and  (2)  the  GRI  2.11  does  not 
contain  the  N2H2,  N2H3,  and  N2H4  species  that  are  present  in  the  SSLA  mechanism.  All  the 
calculations  were  performed  with  both  the  thermal  diffusion  and  the  multicomponent  transport 
package  option,  and  with  the  normal  boundary  conditions  (i.e.,  no  recombination  of  H  atoms  to  Hj 
via  burner  surface  reaction).  Previous  work  under  similar  conditions  showed  that  the  inclusion  of 
this  recombination  effect  only  results  in  a  minor  change  to  the  H  and  H2  profiles  very  close  to  the 
burner  surface.  The  SNL  transport  and  thermochemical  databases  were  used  for  the  SSLA 
calculations,  whereas  the  supplied  transport  and  fliermodynamic  database  default  files  were  used  for 
the  GRI  2.11  calculations.  A  pansal  of  the  thermodynamic  databases  revealed  that  the  (AH°f)298 
values  for  the  same  species  were  within  a  few  tenths  of  kilocalories/mole  of  each  other,  except  for 
NNH.  The  (AH°f)298  values  for  this  species  are  59.6  and  58.5  kcal/mol  in  the  GRI  2.11  and  SNL 
databases,  respectively.  The  effects  of  the  different  databases  on  our  calculations  were  checked  by 
replacing  the  GRI 2.11  recommended  ones  with  those  provided  by  SNL.  The  neat  flame  calculations 
showed  approximately  2  and  1%  increase  in  the  respective  H2  and  N2O  concentrations  at  2.0  mm, 
and  a  0.6  and  0.3  %  decrease  in  the  respective  HjO  and  N2  concentrations  at  15  mm.  For  the  doped 
flame,  the  percent  change  in  concentration  of  the  reactants  is  halved,  whereas  that  of  the  products 
is  doubled  compared  to  the  neat  flame.  The  concentration  of  OH,  NH,  and  0-atom  decreases  from 
4-12%,  with  OH  showing  the  least  change,  for  both  neat  and  doped  flames.  Li  the  neat  and  doped 
flame,  NO  showed  approximately  a  7  and  1 1%  increase  in  concentration,  respectively,  whereas  O2 


7 


showed  a  1 5%  and  2%  decrease  in  concentration.  As  expected,  the  NNH  species  showed  the  largest 
change  in  concentration,  approximately  a  41  %  increase  for  the  neat  flame  and  a  35%  increase  for  the 
doped  flame. 

4.  Results  and  Discussion 

The  modeling  results  of  the  neat  and  NKj-doped  Hj/NjO/Ar  flames  using  the  SSLA  and 
GRI  2.11  mechanisms  reveal  that  the  major  species  profiles  agree  well  with  those  observed 
experimentally,  but  that  the  NH3,  O2,  and  NH  profiles  do  not.  First,  the  modeling  results 
underpredict  the  O2  absolute  concentration  in  the  neat  flame  and  do  not  predict  the  drop  in  O2 
concentration  with  the  addition  of  NH3.  Second,  they  overpredict  NH3  concentrations  throughout 
the  flame  for  the  doped  flame.  Last,  they  do  not  predict  the  decay  of  the  NH  profiles  for  both  neat 
and  NHj-doped  flames.  The  details  in  the  modeled  and  experimental  profiles  is  discussed  in  the 
latter  part  of  this  paper.  Insightful  information  on  which  reactions  to  alter  in  the  mechanisms  for 
better  prediction  of  the  experimental  profiles  is  obtained  by  sensitivity  analyses. 

Presented  in  Table  1  are  the  normalized  sensitivity  coefficients  for  O2  at  15.0  mm  above  the 
burner  surface  for  both  neat  and  NH3-doped  flames  using  both  the  SSLA  and  GRI 2.11  mechanisms. 
The  PREMIX  code  calculates  raw  sensitivity  coefficients  for  each  reaction  and  species  that  are  then 
normalized  according  to  the  following  equation: 

Sik  =  (A/Xi,J(aX^aAi),  (4) 

where  is  the  normalized  sensitivity  coefficient,  Aj  is  the  Arrhenius  A  coefficient  or  reaction  I,  and 
Xif  n,  is  the  maxiiTnim  mole  fraction  of  species  k.  Positive  values  indicate  that  an  increase  in  the 
reaction  rate  results  in  an  increase  in  the  species  mole  fraction,  whereas  negative  values  indicate  that 
an  increase  in  the  reaction  rate  results  in  a  decrease  in  the  species  mole  fraction.  As  revealed  in 
Table  1 , 02  shows  a  very  strong  positive  sensitivity  to  the  N2O  +M =N2 +0  +M  reaction  and  a  negative 
sensitivity  to  the  NH+NO  =N20  +H  reaction  for  both  mechanisms.  O2  also  shows  a  strong  sensitivity 
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Table  1.  Sensitivity  Coefficients  for  O2  in  a  Neat  and  NH3-Doped  Flame  at  15.0  mm  Above  the 
Burner  Surface 


Sensitivity  Coefficients  (Relative)^  || 

SSLA 

_ ^ _ 1 

Sign 

Neat 

Doped 

Neat 

Doped 

+ 

100.0 

100.0 

100.0 

100 

81.6 

16.7 

81.7 

109 

- 

NP 

NP 

44.0 

310 

- 

49.8 

28.0 

22.1 

74 

+ 

33.4 

26.2 

20.7 

98 

- 

30.8 

21.8 

18.7 

75(RR) 

+ 

24.3 

22.9 

18.5 

44 

- 

17.0 

7.5 

9.2 

36 

+ 

16.3 

-101.4 

+37.9 

+173 

+ 

13.9 

<4.5 

14.4 

59 

- 

11.0 

59.1 

19.9 

638 

8.4 

-26.1 

+4.3 

-800 

8.5 

<4.5 

6.8 

2 

- 

4.9 

30.4 

8.7 

740 

+ 

4.3 

17.8 

4.6 

408 

7.5 

16.9 

29 

+ 

4.7 

8.7 

-191 

- 

11.5 

1.3 

154 

— 

<1 

13.8 

<1 

159 

+ 

<1 

8.8 

NP 

NP 

+ 

<1 

6.8 

<1 

-64 

+ 

<1 

4.8 

NP 

NP 

— 

<1 

<1 

+4.7 

+239 

- 

1.5 

1.0 

+2.3 

+77 

— 

1.2 

<1 

-5.6 

<10 

+ 

<1 

<1 

1.9 

155 

— 

<1 

<1 

1.1 

90 

- 

<1 

1.2 

<1 

65 

<1 

1.3 

<1 

18 

+ 

<1 

<1 

<1 

36 

+ 

<1 

<1 

<1 

34 

Reaction 


|N20+M=N2+0+M 

NH+N0=N20+H 

NH+H20=HN0+H2 

NH+OH=HNO+H 

jN0+N=N2+0 

NO+H=N+OH 

NNH+0=N0+NH 

NH+0=N0+H 

N20+H=0H+N2 

NH+N0=N2+0H 

H2+0H=H20+H 

1H+02=0H+0 

N20+0=N0+N0 

H2+0=H+0H 

N20+0=N2+02 

HN0+H=N0+H2 

NO+H+M=HNO+M 

NH2+0=HN0+H 

0H+0H=H20+0 

NH2+N0=NNH+0H 

NH3+0H=NH2+H20 

NH2+NH=N2H2+H 

NH+H=N+H2 

NH+0H=N+H20 

HN0+0H=N0+H20 

N20+0H=N2+H02 

H+H02=0H+0H 

NH3+0=NH2+0H 

NH+02=HN0+0 

H+H02=02+H2 

0H+H02=02+H20 


Notes:  SSLA  =  neat  (100  units  =  0.0454),  doped  (100  units  =  0.323). 

GRI  =  neat  (100  units  =  0.594),  doped  (100  units  =  0.0261). 

NP  =  not  jwesent  in  die  mechanism. 

RR  =  reverse  reaction. 

“  The  sensitivity  coefficients  is  normalized  logarithmically  using  the  maximum  O2  mole  fraction  and  then  scaled  to 
the  N20+M=N2+0+M  reaction. 
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to  the  NH+OH=HNO+H,  N0+N=N2+0,  and  NO+H=N+OH  reactions  for  both  mechanisms.  The 
reaction  NH+H20=HN0+H2,  which  has  an  appreciable  negative  sensitivity  coefficient  for  O2  in  the 
GRI  2.11  mechanism,  is  absent  from  the  SSLA  mechanism.  When  this  reaction  is  included  in  the 
SSLA  mechanism,  the  O2  and  NO  concentrations  decreased  by  56  and  30%  in  the  neat  and  doped 
flame,  respectively,  at  15  mm,  increasing  the  discrepancy  between  modeled  and  observed  values. 

Presented  in  Figure  3  are  the  of  NH3  normalized  sensitivities  as  a  function  of  height  above  burner 
surface  for  the  doped  flame  using  both  the  SSLA  and  GRI  2.11  mechanisms.  NH3  is  very  sensitive 
to  N20+M=N2+0+M,  N20+H=0H+N2,  H2+0H=H20+H,  and  NH3+0H=NH2+H20  for  bothSSLA 
and  GRI  2.11  mechanisms.  NH3  also  shows  a  negative  sensitivity  to  reactions  NH2+NO=NNH+OH 
and  NH2+NH=N2H2+H  and  apositive  sensitivity  to  OH+NH=HNO+H  throughout  the  flame  for  the 
SSLA  mechanism.  In  contrast,  NH3  is  not  sensitive  to  the  OH+NH=HNO+H  reaction  for  the 
GRI  2.11  mechanism.  The  NH3  sensitivity  for  reactions  NH2+NH=N2H2+H  and 
NH2+N0=NNH+0H  is  not  displayed  for  the  GRI 2.11  mechanism  because  they  are  absent  from  the 
mechanism.  The  NH2+N0=N2+H20  reaction  is  also  not  included  in  the  GRI  2.11  mechanism. 
Another  difference  between  the  two  mechanisms  shown  in  Figure  3  is  that  all  the  NH3  sensitivity 
coefficients  for  the  SSLA  mechanism  peak  at  approximately  5  mm  and  then  gradually  decrease, 
reaching  ~'10%  of  their  values  near  25  mm  In  contrast,  NH3  shows  substantial  sensitivity  to  many 
reactions  above  25  mm  for  the  GRI  2.11  mechanism  Little,  if  any,  NH3  is  observed  experimentally 
above  10  mm 

The  O2  and  NH3  species  are  very  sensitive  to  the  N2O  +M=N2 +0 +M  reaction  as  shown  in  T able  1 
and  Figure  3,  respectively,  and  an  increase  in  its  rate  expression  will  increase  the  postflame 
concentration  of  O2  and  decrease  that  of  NH3.  This  trend  is  in  agreement  with  the  experimental 
results.  An  increase  in  the  N20+M=N2+0+M  reaction’s  rate  expression  also  decreases  the  NH 
concentration  for  the  neat  and  NH3-doped  flames  near  the  burner  smface,  as  shown  in  Table  2,  and 
better  predicts  its  decay  in  both  flames.  O2  is  also  very  sensitive  to  the  NH+N0=N20+H  reaction. 
A  decrease  in  its  rate  constant  results  in  an  increase  in  the  O2  postflame  concentration,  consistent  with 
the  e3q)erimental  results.  The  NH3  concentration  throughout  the  flame  would  not  change,  however, 
because  NH3  is  insensitive  to  this  reaction.  Table  1  also  shows  that  O2  is  sensitive  to 
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nh3  sensitivity  coefficients 
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Figure  3.  Plot  of  NH3  Sensitivity  Coefficients  as  a  Function  of  Hdght  Above  Burner  for 
Various  Reactions  Using  Both  SSLA  (Top)  and  GRI  2.11  (Bottom)  Mechanisms. 
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Table  2.  Sensitivity  Coefficients  for  NH  in  a  Neat  and  NHj-Doped  Flame  at  4.1  and  53  mm, 
Respectively,  Above  the  Burner  Surface 


Reaction 

Sensitivity  Coefficients  (Relative)®  || 

SSLA 

GRI  1 

Sign 

Neat 

Doped 

Neat 

Doped 

NH+N0=N20+H 

+ 

100.0 

99.9 

84.5 

55.7 

N20+H=0H+N2 

92.5 

107.1 

100.0 

100.0 

NH+H20=HN0+H2 

- 

NP 

NP 

28.5 

57.0 

N20+M=N2+0+M 

- 

37.8 

14.8 

23.6 

8.2 

N+H2=NH+H 

- 

33.3 

64.1 

22.9 

29.4(RR) 

NH+OH=HNO+H 

- 

27.6 

45.6 

16.1 

17.9 

H2+0H=H20+H 

+ 

23.3 

22.9 

19.3 

3.6 

NH+N0=N2+0H 

- 

16.1 

18.5 

14.9 

16.1 

NNH+0=NO+NH 

- 

14.2 

12.0 

8.2 

5.4 

NH+0H=N+H20 

12.7 

18.2 

11.5 

12.0 

NH+NH=N2+H2 

- 

9.2 

42.3 

NP 

NP 

NH+0=NO+H 

- 

8.2 

11.3 

5.0 

8.0 

NO+H=N+OH 

- 

4.4 

11.4 

3.5 

8.6(RR) 

N0+N=N2+0 

+ 

4.3 

18.4 

3.1 

8.0 

H2+0=H+0H 

+ 

2.8 

11.9 

2.0 

10.1 

NH2+H=NH+H2 

+ 

2.7 

30.5 

2.0 

23.8 

NH2+NH=N2H2+H 

- 

<1 

41.5 

NP 

NP 

NH2+0H=NH+H20 

+ 

<1 

21.0 

<1 

12.8 

NH3+0H=NH2+H20 

+ 

<1 

16.2 

<1 

19.6 

NH2+0=HN0+H 

- 

<1 

10.3 

<1 

7.4 

NH2+N=N2+H+H 

- 

<1 

7.8 

NP 

NP 

NH3+H=NH2+H2 

- 

<1 

<1 

<1 

9.8 

HN0+H=H2+N0 

— 

— 

— 

-4.5 

-24.7 

H+NO+M=HNO+M 

— 

— . 

— 

-5.2 

-15.7 

Notes:  NP  =  not  present  in  flie  mechanism. 

RR  =  reverse  reaction. 

*  The  sensitivity  coefficients  are  normalized  logarithmically  using  the  maximum  NH  mole  fraction  and  then  scaled  to 
the  NH+N0=N20+H  reaction  for  the  SSLA  mechanism  (neat  [100  units  =  0.518],  doped  [100  units  =  0.0341])  and 
N20+H=0H+N2  reaction  for  the  GRI  mechanism  (neat  [100  units  =  -0.556],  doped  [100  units  =  -0.365]). 


reactions  NH+OH=HNO+H,  NO+N  =N2+0,  and  NO+H=N+OH,  and  their  rate  coefficients  could 
be  altered  to  increase  the  Oj  postflame  concentration.  However,  these  reactions  affect  the  Oj 
postflame  concentration  in  NHj/NjO  flame  in  a  similar  manner,  and  the  desired  result  would  be 
contrary  to  what  is  observed  experimentally  [17].  Thus,  the  rate  expressions  for  these  reactions  are 
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not  altered.  NH3  is  also  sensitive  to  the  N20+H=0H+N2  reaction  (see  Figure  3).  A  decrease  in  its 
rate  constant  produces  a  lower  NH3  concentration  throughout  the  flame,  in  agreement  with  the 
experimental  results,  but  also  a  lower  O2  postflame  concentration,  contrary  to  what  is  observed 
experimentally.  Thus,  the  rate  constant  for  this  reaction  is  also  not  altered. 

The  SSLA  mechanism  was  modified  by  altering  the  rate  expression  and  third  body  H2O  efficiency 
for  the  N20+M=N2+0+M  reaction  and  the  rate  expression  for  the  NH+N0=N20+H  reaction.  The 
low-pressure  limit  rate  expression  of  the  N2O +M=N2 +0 +M  reaction  for  N2  in  the  SSLA  mechanism 
is  that  reported  in  Rohrig  et  al.  [18],  k  =  5.97  x  10*^  exp(-56640/T),  with  a  third  body  H2O/N2 
efficiency  ratio  of  7.5  [19].  In  contrast,  the  rate  expression  in  the  GR/ 2.2/  mechanism  for  Nj  is  that 
reported  by  Glarborg  and  coworkers,  k  =  6.2  x  10^“*  exp(-  56100/T)  with  a  H2O/N2  ratio  of  6.0  [20]. 
This  expression  is  in  excellent  agreement  with  that  reported  by  Hanson  and  Salimian  [21].  At  the 
peak  ten5)erature  in  the  H2/N20/Ar  flame,  ~2000  K,  the  GRI2JI  value  is  approximately  26%  larger 
than  that  used  in  the  SSLA  mechanism.  The  H2O  and  N2  third-body  efficiency  values  in  the  SSLA 
mechanism  are  those  reported  by  Glarborg  and  co workers  [19].  Their  values  relative  to  Ar  are 
12  ±3.5  and  1.7  ±0.3,  respectively,  in  the  1,000-1,400-K  range.  The  SSLA  mechanism  is  thus 
modified  to  include  the  rate  constant  expression  from  the  GRI  2.11  mechanism,  with  an  effective 
H2O/N2  ratio  of  1 1 ,  the  upper  value  of  experimental  uncertainty  reported  by  Glarborg  and  co  workers 
[19].  A  similar  effect  in  the  modeling  results  can  be  achieved  by  using  the  rate  constant  reported  in 
Rohrig  et  al.  [18]  with  a  larger  effective  ratio  for  H2O/N2,  which  is  certainly  plausible  for  our  flame 
conditions. 

The  rate  expression  for  the  NH+NO=N20+H  reaction  in  the  SSLA  mechanism  is  that  obtained 
by  Miner  and  Melius  [22]  from  BAC-MP4  calculations,  k  =  2.94  x  10^^  T'®'”  -  2.16  x 
In  the  SSLA  modified  mechanism,  it  is  decreased  by  25%,  certainly  within  the  theoretical  uncertainty. 
For  con^arison  purposes,  the  rate  expression  in  the  GRI  2.11  mechanism,  k  =  4.16  x  10^“*  j-0.45^ 
yields  a  rate  constant  at  2,000  K,  which  is  approximately  32%  greater  than  that  used  in  the  SSLA. 
This  reaction  strongly  affects  both  the  SSLA  and  GRI  2.11  calculated  postflame  NO  concentration. 
Calculations  for  the  neat  flame  using  both  mechanisms  show  that  the  NO  sensitivity  coefficient  for 
this  reaction  is  positive  and  is  second  (55%)  to  that  of  the  N20+H=0H+N2  (100%)  reaction  in 
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absolute  value.  In  the  doped  flame,  the  SSLA  and  GRI  2.11  sensitivity  coefficients  for  the 
NH+N0=N20+H  reaction  are  ~33%  less  than  that  of  the  N20+H=0H+N2  reaction. 

Four  different  sets  of  calculations  are  performed  for  neat  and  NHj-doped  flames:  ( 1 )  SSLA  -  uses 
the  SSLA  mechanism  and  the  measured  tenqjerature  profile  as  input;  (2)  GRI  2.11  -  uses  the 
GRI  2.11  mechanism  and  the  measured  ten:5)erature  profile  as  input;  (3)  SSLA  modified  -  uses  the 
SSLA  mechanism  in  which  the  N20+M=N2+0+M  and  NH+N0=N20+H  reactions  are  altered,  as 
discussed  previously,  and  a  5%  increase  in  the  temperature  profile  as  input  (twice  the  standard 
deviation  of  our  experimental  uncertainty);  and  (4)  GRI  2.11  modified  -  uses  the  GRI  2.11 
mechanism  in  which  the  aforementioned  two  reactions  are  altered  as  in  the  SSLA  mechanism  and  a 
5%  increase  in  the  temperature  profile  as  input.  The  modified  GRI 2.11  mechanism  is  not  optimized, 
and  the  calculations  are  included  for  comparison  purposes. 

The  H2,  N2O,  H2O,  and  N2  calculated  profiles  for  both  the  neat  and  NHj-doped  H2/N20/Ar  flames 
agree  well  with  those  obtained  experimentally  using  all  four  mechanisms.  The  calculated  NH3  profiles 
for  the  doped  flame,  along  with  the  experimental  profile,  are  shown  in  Figure  4.  The  experiment  and 
model  calculations  show  a  negligible  mole  fi-action  for  the  NH3  in  the  neat  flame.  The  SSLA 
modified  model  overpredicts  the  NH3  by  6%  near  the  burner  surface,  which  is  within  the  experimental 
uncertainty  of  the  measurement.  The  model  also  indicates  that,  at  the  height  of  15  mm  (postflame 
region),  the  NH3  dopant  decreases  to  a  level  of  0.6%  of  the  added  amount,  similar  to  that  observed 
experimentally.  In  contrast,  the  SSLA  calculations  reveal  that  the  calculated  NH3  mole  fi'action  (at 
1 5  mm)  decreases  to  3%  of  the  added  amount.  Calculations  using  the  modified  GRI  mechanism  also 
predict  this  result.  The  NH3  profile  calculated  with  the  GRI  mechanism  has  the  maximum  deviation 
firom  the  experimental  results,  as  shown  in  Figure  3,  and  shows  15%  of  NH3  remaining  unconsumed 
in  the  postflame  region. 

Rate  analysis  for  the  SSLA  modified  calculations  reveal  that  all  of  the  NH3  is  consumed  by  the 
NH3+0H=NH2+H20,  NH3+H=NH2+H2,  and  NH3+0=NH2+0H  reactions,  the  first  accounting  for 
~86%  of  its  consunq)tion.  The  reactions  NH3+M=NH2+H+M,  NH2+NH2=NH+NH3,  and 
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figure  4.  Experimental  and  Modeled  NH3  Profiles. 

HN0+NH3=NH3+N0  account  for  most  of  its  formation,  but  their  total  rate  is  ~3.8%  of  the 
consumption  rate. 

The  experimental  and  modeled  O2  profiles  for  the  neat  and  NH3-doped  flames  are  shown  in 
Figure  5.  The  MB/MS  measurements  reveal  a  90%  decrease  in  Oj  in  the  burnt  gas  region  of  the 
NHj-doped  flame.  The  SSLA  modified  modeling  results  show  that  the  amount  of  O2  formed  in  the 
doped  flame  decreases  by  86%  and  underpredicts  the  O2  formation  in  the  neat  flame  by  5%  at 
15  mm.  The  modified  GRI  2.11  mechanism  also  underpredicts  the  amount  of  O2  in  the  neat  flame 
by  13%,  but  shows  a  relative  drop  in  the  doped  flame  of  90%.  Overall,  both  of  the  aforementioned 
mechanisms  well  predict  the  shape  of  the  O2  profiles.  However,  the  modified  GRI  2.11  results 
display  a  more  pronounced  decay  in  O2  profile  in  the  burnt  gas  region  of  the  doped  flame.  The  SSLA 
and  GRI  2.11  display  the  maximum  deviation  (39%  and  57%,  respectively)  from  the  measured  O2 
mole  fraction  in  the  neat  flame. 
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Figure  5.  Experimental  and  Modeled  Oj  Profiles. 

I*resented  in  Figure  6  are  the  major  reactions  forming  and  consuming  Oj,  along  with  the  total  O2 
rate,  for  both  neat  and  NHj-doped  flames.  For  the  neat  flame,  all  of  the  O2  is  produced  by  the 
N20+0=N2+02  and  H+02=0H+0  reactions.  Reactions  N0+0=N+02  and  NH+02=HN0+0 
contribute  slightly  toward  its  consumption,  whereas  the  02+H+M=H02+M  reaction  has  little  effect 
on  both  its  formation  or  consumption.  For  the  doped  flame,  the  N20+0=N2+02  reaction  also  forms 
O2  and  die  H+02=0H+0  reaction  forms  O2  up  to  7.5  mm  and  consumes  it  at  distances  greater  than 
7.5  mm.  Reactions  NO  +0  =N+02  and  NH  +O2 =HNO  +0  contribute  slighdy  toward  O2  consumption, 
as  in  the  neat  flame,  but  are  approximately  50%  less  effective.  The  net  integrated  rate  of  O2 
production  (0-20  mm)  in  the  neat  flame  is  thus  ~10  times  larger  than  that  in  the  doped  flame 
primarily  because  the  N2O  +0=N2+02  reaction  has  a  smaller  positive  contribution  in  the  doped  flame 
and  because  the  H +O2  =OH +0  reaction  has  a  smaller  positive  contribution  at  distances  <7.5  mm  and 
a  negative  contribution  at  distances  >7.5  mm. 
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Figure  6.  Oj  Reaction  Rates  From  SSLA  Modified  Calculations  for  Both  Neat  and  NHj-Doped 
Hj^jO/Ar  Flames. 
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It  should  be  noted  that  the  rate  expression  for  the  N20+0=N2+02  reaction  used  in  the  SSLA 
modified  mechanism  is  that  reported  by  Hanson  and  coworkers  [23],  k  =  1.40  x  10'^ 
exp(-  10,800/T).  Altering  this  rate  expression  to  that  recommended  by  Hanson  and  Salimian  [21], 
k  =  1.0  X  lO*'*  exp(-28,020/T)  results  in  a  negligible  change  in  the  computed  Oj  and  NO  mole 
fractions  for  both  neat  and  NH3-doped  flames  at  15  mm.  In  addition,  the  net  integrated  rate  of  O2 
production  for  both  neat  and  doped  flames  is  hardly  changed.  The  peak  rate  for  this  reaction  is, 
however,  decreased  by  approximately  32  and  25%  for  neat  and  doped  flames,  respectively,  and  it 
shifted  approximately  1  mm  away  from  the  burner  surface  for  both  flames. 


Figure  7  shows  the  NO/Ar  ratio  measured  by  MB/MS  in  the  neat  and  doped  flame.  For 
comparison  purposes,  the  NO  profiles  obtained  with  UDF  are  converted  to  NO/ Ar  ratio  using  the 
measured  temperatures  and  the  calculated  Ar  mole  fractions.  The  NO/Ar  (LJF)  profile  is  then 
normalized  to  the  calibrated  MB/MS  value  in  the  neat  flame.  The  UF  from  NO  is  also  quantified 
by  calibrating  the  NO  fluorescence  signals  using  LDF  measurements  from  calibrated  mixtures  of  NO 
and  Ar  at  room  temperature.  The  analysis  incorporates  the  change  in  the  Boltzmann  fraction  of  the 
ground  level  due  to  the  temperature  difference  between  the  calibration  and  die  flame  conditions.  The 
T  .TF  value  is  approximately  50%  lower  than  the  MB/MS  value,  within  the  experimental  uncertainty 
of  the  UF  measurement.  A  similar  NO  UF  uncertainty  is  reported  by  Heard  et  al.  [24]  for  a  30-T orr 
CH4/Air  flame.  The  MB/MS  measurements  show  a  decrease  of  32%  in  the  NO  when  NH3  is  added. 


whereas  the  UF  measurements  show  a  decrease  of  35%.  Considering  the  experimental  uncertainties 
involved,  both  the  UF  and  MB/MS  measurements  show  the  same  level  of  NO  depletion  in  the  doped 
flame.  The  SSLA  modified  calculations  predict  a  38%  drop  in  the  level  of  NO  between  the  neat  and 
NH3-doped  flames,  in  excellent  agreement  with  that  observed  experimentally.  The  calculations  also 
underpredict  slightly  the  NO  mole  fraction  in  both  the  neat  and  doped  flame.  All  of  the  calculations 
used  in  this  study  well  predict  the  shape  of  the  NO  profiles  and  the  mole  fraction  of  NO  in  the  neat 
flame.  The  relative  decrease  in  NO  with  NH3  doping  is  well  predicted,  however,  only  by  the  SSLA 


modified  and  GRI  2.11  models. 


Figure  8  shows  key  reactions  contributing  to  NO  production  and  consumption,  along  with  the 
total  rate,  for  both  neat  and  NH3-doped  flames.  As  expected,  the  peak  rate  of  production  at 
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Figure?.  Experimental  and  Modeled  NO  Profiles  for  Neat  and  NHj-Doped HjWjO/Ar 
Flames.  The  Experimental  Profiles  Are  Obtained  by  MB/MS  (•)  and  LIF  (o). 

approximately  5  mm  corresponds  to  die  maximum  temperature  rise.  The  net  integrated  rate 
(0-30  mm)  of  NO  production  for  the  neat  flame  is  approximately  2.4  x  10"°®  mol/cm^-  s.  Reactions 
NH+N0=N20+H,  HN0+H=N0+H2,  and  NO+H=N+OH  account  for  86%  of  the  NO  production 
flux,  whereas  reactions  NO+N =N2+0,  NH+NO=N2+OH,  and  NNH+0=NO+NH  account  for  over 
99%  of  die  NO  consumption  flux.  The  addition  of  NHj  to  the  neat  flame  results  in  a  net  integrated 
rate  of  1.46  x  10'°®  mol/cm^-s,  a  decrease  of  ~42%  from  the  neat  flame.  This  decrease  is  primarily 
from  the  NH+N0=N20+H  reaction.  The  HNO+H=NO+H2  and  NO+H=N +OH  reactions  produce 
more  NO  in  the  doped  flame  compared  to  the  neat  flame.  Consumption  of  NO  from  reactions 
N0+N=N2+0,  NH+N0=N2+0H,  and  NNH+0=N0+NH  is  also  greater  in  the  doped  flame 
compared  to  die  neat  flame.  In  addition,  the  integrated  reaction  rates  for  the  NH2+NO=NNH+OH 
and  NH2+N0=N2+H20  reactions  are  approximately  20  times  greater  in  the  doped  flame  compared 
to  the  neat  flame,  accounting  for  ~16%  of  the  consumption  of  NO  in  the  doped  flame. 


19 


RATE  (mol/cm^-sec) 


Rt:  Total  Rate  R6:  NNH+0=N0+NH 

R1:  NH+N0=N20+H  R7:  NH+0=N0+H 

R2:  N0+N=N2+0  R8:  N20+0=N0+N0 

R3:  NO+H=N+OH  R9:  NO+H+M=HNO+M 

R4:  HN0+H=N0+H2  R10:  NH2+NO=NNH+OH 

R5:  NH+N0=N2+0H  R1  1 ;  NH2+NO=N2+H 
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Figure  8.  Reaction  Rates  of  NO  From  SSLA  Modified  Calculations  for  Both  Neat  and 
NH3-Doped  H^/NjO/Ar  Flames. 
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Figure  9  shows  the  calculated  and  the  experimental  NH  profiles.  The  experimental  and 
calculated  relative  NH  profiles  are  normalized  to  unity  for  the  neat  flame.  The  modified  SSLA 
model  predicts  a  126%  increase  in  NH  when  NH3  is  added  to  the  flame.  This  agrees  well  with  the 
109%  increase  measured  by  LIF.  The  SSLA  model  calculations  pr^ict  an  89%  increase  in  NH  with 
NH3  doping,  whereas  the  GRI  2.11  and  modified  GRl  2.11  calculations  predict  a  95%  and  159% 
increase,  respectively.  The  predicted  NH  peaks  are  at  3.75  mm  for  the  neat  and  4.75  mm  for  the 
doped  flames,  in  agreement  with  the  experimental  peaks  at  3.5  and  4  mm  for  the  neat  and  doped 
flame,  respectively.  The  SSLA  modified  calculations  well  predict  die  shape  of  NH  profile  in  the  neat 
flame,  including  the  position  of  its  peak  concentration.  The  calculations  also  well  predict  the  relative 
increase  in  NH  concentration  as  the  neat  flame  is  doped  with  NH3  For  the  doped  flame,  all  the 
calculations  show  a  broader  NH  peak  than  that  observed  experimentally,  with  the  SSLA  modified 
results  agreeing  best. 


Figure  9.  NH  Experimental  and  Modeled  Profiles  for  Neat  and  NHj-Doped 

Flames.  The  Experimental  Concentration  Profiles  Are  Converted  to  Mole  Fraction 
Using  the  Experimental  Temperature  Profile.  All  the  Profiles  for  the  Neat  Flame 
Are  Normalized  to  Unity.  Their  Peak  Mole  Fraction  Values  Are  2251  x  10'^ 
1.972  X  10-«,  1.816  X  10 and  1.393  x  10-^  for  the  SSLA,  GRl  2.11,  SSLA  Modified, 
and  GRI  2.11  Modified  Calculations,  Respectively. 
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The  SSLA  modified  calculations  rate  analyses  results  are  presented  in  Figure  10.  The  integrated 
rates  (0-30  mm)  for  the  neat  flame  show  that  almost  all  of  the  NH  is  formed  from  the 
NH+N0=N20+H  reaction  and  64%  is  consumed  from  tiie  N+H2=NH+H,  NH+OH=HNO+H,  and 
NH+N0=N2+0H  reactions.  For  the  doped  flame,  the  NH+N0=N20+H  reaction  still  plays  a  key 
role  in  NH  formation.  However,  the  NH2+H=NH+H2  and  NH2+0H=NH+H20  reactions  are  also 
important  and  account  for  >-48%  of  the  NH  production.  In  contrast,  these  reactions  form  <1%  NH 
in  the  neat  flame.  The  reactions  N+H2=NH+H,  NH+OH=HNO+H,  NH+NO=N2+OH  and 
NH+OH=N +H2O  are  important  for  NH  consumption  in  the  doped  flame,  as  in  the  neat  flame,  with 
the  integrated  (0-30  mm)  N  +H2=NH+H  reaction  rate  being  a  factor  of  2  larger  than  that  of  the  neat 
flame. 

The  experimental  and  modeled  OH  profiles  for  the  neat  and  NHj-doped  flame  are  shown  in 
Figure  11.  All  the  calculated  OH  profiles  for  the  neat  flame  were  also  normalized  to  unity.  The 
shape  of  the  OH  profiles  are  predicted  very  well  by  the  SSLA  modified  model.  The  calculations 
show  a  53%  drop  in  the  OH  in  the  burnt  gas  region  of  the  flame  when  NH3  is  added,  which  agrees 
rather  well  with  the  experimental  decrease  of  43%.  As  shown  in  Figure  11,  using  the  other  three 
mechanisms  results  in  the  OH  profile  decaying  in  postflame  region,  with  the  GRl  2.11  results  being 
the  most  pronounced.  These  mechanisms  also  overpredict  the  percent  decrease  in  the  amount  of  OH 
with  NH3  doping.  Rate  analyses  of  the  SSLA  modified  calculations  reveal  that  in  the  neat  flame  OH 
is  formed  primarily  from  the  N20+H=0H+N2  reaction  and  consumed  by  the  H2+0H=H20+H 
reaction,  as  expected.  The  H2+0=H+OH  reaction  also  plays  a  lesser,  but  an  important  role  in  die 
OH  formation.  For  the  doped  flame,  the  aforementioned  reactions  also  play  key  roles  in  the 
production  and  consumption  of  OH.  The  NH3+0H=NH2+H20  also  accoimts  for  '-9%  of  the  OH 
consumption  in  the  doped  flame,  compared  to  the  neat  flame  in  which  it  plays  a  minor  role. 

The  experimental  and  calculated  relative  concentration  profiles  for  the  0-atom  are  shown  in 
Figure  12.  Overall,  the  models  fairly  well  predict  the  shape  of  the  O-atom  profile  in  the  neat  flame. 
However,  none  of  the  models  predict  the  10%  decay  observed  in  die  O-atom  concentration  profile 
in  the  postflame  region.  In  fact,  all  of  the  calculated  profiles  exhibit  a  plateau  in  this  region.  The 
modified  SSLA  calculations  show  the  best  agreement  in  the  decrease  in  the  O-atom  concentration 
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Figure  10.  NH  Reaction  Rates  From  SSLA  Modified  Calculations  for  Both  Neat  and 
NHj’Doped  H/NjO/Ar  Flames. 
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Figure  11.  OH  Experimental  and  Modeled  Profiles  for  Neat  and  NH3-Doped  Hj/NjO/Ar 
Flames.  The  Experimental  Concentration  Profiles  Are  Converted  to  Mole  Fraction 
Using  the  Experimental  Temperature  Profile.  All  the  Profiles  for  the  Neat  Flame 
Are  Normalized  to  Unity.  Their  Mole  Fraction  Values  at  15  mm  Are  9.932  x  10"^, 
4.140  X 10 7.211  X 10®,  and  6.689  x  10®  for  the  SSLA,  GRI2J1,  SSLA  Modified, 
and  GRI  2.11  Modified  Calculations,  Respectively. 

in  the  postflame  region  when  NH3  is  added.  A.  11%  decrease  is  calculated  compared  to  a  68% 
decrease  measured  by  LIF.  The  models  do  not  very  well  predict  the  shape  of  the  0-atom  profile  in 
the  doped  flame.  All  of  the  mechanisms  used  in  the  modeling  of  the  flames  show  a  slower  chemistry 
than  experimentally  observed  (i.e.,  the  calculated  profiles  are  shifted  away  from  the  burner) 
especially  in  the  doped  flame.  One  possible  explanation  of  this  effect  is  photochemical  interference 
from  0-atom  [25].  Our  experimental  findings  do  not  support  this  explanation,  however,  because 
there  should  be  a  decrease  and  not  an  increase  in  the  interference  effects  in  the  doped  flame,  where 
the  amount  of  O2  present  is  greatly  reduced.  Alternative  explanations  include  photochemical  effects 
for  other  O-atom-containing  species  and/or  further  refinements  of  the  modified  SSLA  mechanism. 
Experiments  and  PREMIX  calculations  along  these  lines  are  ongoing. 
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Figure  12.  0-Atom  Experimental  and  Modeled  Profiles  for  Neat  and  NHj-Doped 

Flames.  The  Experimental  Concentration  Profiles  Are  Converted  to  Mole  Fraction 
Using  the  Experimental  Temperature  Profile.  All  the  Profiles  for  the  Neat  Flame 
Are  Normalized  to  Unity.  Their  Mole  Fraction  Values  at  15  mm  Are  5.125  x  lO''', 
3.826  X 10 9.180  x  10 and  SJ2411  x  lO  "*  for  the  SSLA,  GRI2J1,  SSLA  Modified, 
and  GRI  2.11  Modified  Calculations,  Respectively. 

5.  Conclusion 

A  combined  experimental  and  modeling  study  of  neat  and  NHj-dqped  Hj/NjO  flames  has  been 
performed.  Species  concentrations  were  measured  by  MB/MS  and/or  IIF,  and  flame  temperatures 
were  measured  by  bofli  OH  and  NH  LIF  and  thin- wire  thermometiy .  The  experimental  profiles  were 
compared  widi  calculated  profiles  generated  using  PREMIX  with  the  following  detailed  chemical 
mechanisms:  (1)  SSLA,  (2)  modified  SSLA,  (3)  GRI  2.11,  and  (4)  modified  GRI  2.11.  Overall,  the 
major  experimental  species  profiles  for  both  neat  and  NHj-doped  flames  are  predicted  fairly  well 
using  all  four  mechanisms.  The  minor  species  profiles  are  modeled  best,  however,  using  the 
modified  SSLA  mechanism.  In  this  mechanism  the  rate  expression  for  the  NH+N0=N20+H 


25 


reaction  is  decreased  by  25%,  the  rate  expression  for  the  N20+M=N2+0+M  reaction  is  increased 
25%,  and  the  H2O  third-body  efficiency  for  the  latter  reaction  increased  ~47%,  the  upper  limit  of  the 
experimental  uncertainty.  The  modified  SSLA.  calculations  rather  well  predict  the  species  profiles 
for  both  neat  and  doped  flames,  as  well  as  the  increase  in  postflame  O2,  NO,  OH,  and  0-atom 
concentrations,  and  the  increase  in  the  peak  NH  concentration  with  the  addition  of  4%  NH3.  The 
calculations  show  a  decrease  of  38%  in  NO,  36%  in  O2,  53%  in  OH,  and  77%  in  the  0-atom,  and 
an  increase  of  126%  in  NH,  compared  to  a  decrease  of  32  and  35%  in  NO  measured  by  UF  and 
MB/MS,  respectively,  and  a  decrease  of  90%  in  O2,  43%  in  OH,  and  68%  in  the  0-atom  and  an 
increase  of  109%  in  NH  measured  by  LIF.  The  NH  and  0-atom  profiles  for  the  neat  flame  are  also 
predicted  the  best  and  fairly  well  by  the  modified  SSLA.  For  the  doped  flame,  however,  all 
mftrhanisms  predict  a  broader  NH  profile  and  an  0-atom  profile  that  is  shifted  ~2  mm  away  from 
the  burner  compared  to  that  measured  experimentally.  These  discrepancies  are  presently  under 
investigation. 
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